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Abstract 

We present a dynamic nonlocal hybrid Monte Carlo algorithm consisting of pivot 
and "cut-and-permute" moves. The algorithm is suitable for the study of polymers 
in semiconfined geometries at the ordinary transition, where the pivot algorithm 
exhibits quasi-ergodic problems. The dynamic properties of the proposed algorithm 
are studied in d = 3. The hybrid dynamics is ergodic and exhibits the same optimal 
critical behavior as the pivot algorithm in the bulk. 



1 Introduction 



In this paper we introduce and discuss the properties of a hybrid Monte Carlo algorithm 
which can be used to study the equilibrium properties of a polymer molecule grafted to a 
surface. 

The system can be experimentally obtained either by chemically grafting one polymer 
end in an irreversible way, or, as for surfactants, via physical adsorption of an endgroup or 
one of the two blocks in a diblock copolymer. In this last case, the process is reversible, 
since the attached end can desorb both when the temperature increases and when the 
solvent quality changes. 

Systems of this kind include polymers grafted at colloidal particles or surfaces in so- 
lution which can help stabilizing against flocculation IQ or polymers grafted at flexible 
lipid membranes exerting on the membrane a bending force which is proportional to the 
monomer concentration at the membrane 

In what follows we will consider a polymer in a good solvent which is grafted at a surface 
S and interacts repulsively with S. We will focus in particular on the case in which is a 
fiat surface. The polymer will be modelled by an iV-step lattice self-avoiding walk (SAW), 
which provides a good description of the critical behavior of polymer molecules in the bulk 
as well as in confined geometries 

Efficient simulations of SAWs can be obtained by using nonlocal algorithms. For in- 
stance, the pivot algorithm (see Refs. ||8|-[TT| and references therein) is optimal, up to a 



constant factor, for sampling global observables in the fixed-A^ free-endpoint ensemble in 
the bulk since the autocorrelation time in CPU units is simply proportional to A^. In the 
presence of surfaces, this algorithm is not as efficient and in some cases it is not even er- 
godic (this is the case of a two-dimensional strip, see Ref. |12|). In the presence of a single 
impenetrable plane, the algorithm is ergodic, but still one expects it to be inefficient since 
the initial part of the walk will be rarely updated. For instance, we will show that, in 
order to update the direction of one of the first links pointing along the normal vector to 
the surface, the walk should have an extremely unlikely geometrical shape. This means 
that, in order to explore the relevant phase space, the algorithm has to go through highly 
improbable configurations. Therefore, even if ergodic, the algorithm has an autocorrelation 
time that grows rapidly with A^. 

To overcome these problems, we propose here a hybrid algorithm based on the pivot 
move and on a cut-and-permute move. This last move has been introduced in Ref. in a 
hybrid algorithm working in the fixed-endpoint, variable- A^ ensemble. It consists in cutting 
the walk into two parts and in rebuilding it by concatenating the two parts in the opposite 
order. As we will discuss, such a move is quite efficient in changing the configuration of 
the walk near the grafted endpoint. Thus, if we combine the pivot move and the cut-and- 
permute move, we obtain an algorithm which does not have the quasi-ergodicity problems 
of the pivot algorithm. Also, by a careful implementation, it is possible to obtain the same 
optimal behavior as the pivot algorithm for polymers in the bulk: the autocorrelation times 
of global observables in CPU units increase only as A^. It must be noted that, even though 
we consider here only the case of a fiat surface, the results should apply to any convex 
surface, for instance to polymers grafted outside a sphere: in this respect, the plane should 
be the worst case. Note also that we do not consider here any interaction between the 
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surface and the walk other than the excluded-volume interaction. Nonetheless, the results 
should also hold in the presence of attractive interactions as long as they are sufficiently 
weak and the walk is not absorbed. In this case, it is probably important to add local 
moves (and, perhaps, bilocal moves as defined in Ref. |jl2|) to speed up the dynamics near 
the surface. Finally, we want to notice that this algorithm is also needed if we want to 



apply the join-and-cut algorithm of Ref. [14| to walks in the presence of a surface. Indeed, 



the ideas that are presented here apply directly to that algorithm, so that by using the 
cut-and-permute together with the pivot move, one should be able to have a version of the 
join-and-cut algorithm which works reasonably well also in the presence of a surface. 

The paper is organized as follows. In Sec. 2 and 3 we consider the pivot and the cut- 
and-permute move and discuss in detail the acceptance fraction, paying particular attention 
to those moves in which the pivot or the cutting point is near the grafted endpoint of the 
walk. In Sec. 4 we discuss the implementation of the two moves and compute the scaling 
behavior of the average CPU time for each move. In Sec. 5 we discuss the full algorithm, 
showing its optimal behavior. 



2 The Pivot move in the presence of a surface 

In this paper we will consider iV-step self-avoiding walks on a d-dimensional lattice if" in 
the presence of an excluded surface S of equation S{x) = 0, x G if" . A SAW is a sequence 
of lattice points {cuq, . . . -.ojj^^, such that Wj, tUj+i are lattice neighbours and uJi ^ ojj for 
i ^ j- The walk is confined in the outward half-space S{x) > with its first vertex grafted 
at S and fixed at position ujq. The state space is therefore 

^%{uJo) = {uj = {uo, . . .,ujn}\S (uJi) > , i = l,...,N anduJi^ujj , j}. (1) 

As probability measure on the space Qj^f[uJo), we will consider the uniform one, which 
gives equal weight to every allowed walk, i.e. vr (tu) = where is the cardinality of 

O|r(a;o). In the following we will assume S* to be a (d — 1) hyperplane of equation z = 0. 

In Qff{u!o) we will first consider the pivot algorithm The elementary move of the 

algorithm in the bulk consists of the following steps: 

1. Given a walk configuration u = {ujo, ■ ■ ■ , wn}, an integer A; G {0, . . . , — 1} is chosen 
at random and the corresponding monomer u;^ is taken as pivot point. 

2. An element g of the symmetry group of the lattice is chosen at random with probabil- 
ity Pg and a new walk is built by applying it to the part of the chain which follows the 
pivot point, using uj^ as fixed point of the transformation. A walk u' = {u'q, . . . , cu^}, 
where the new monomer coordinates are given by 

, ^ f for 1 < i < fc, 

' \uJk + 9{uJi -uJk) for A; + 1 < z < AT , ^"^^ 

is obtained. 
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3. If the walk uo' is not self- avoiding the move is rejected and the original walk uj is kept. 
Otherwise, uo' is taken as new current walk. 

We require Pg = Pg-i in order to satisfy detailed balance. In the bulk the algorithm is 
ergodic if all axes reflections and either all 7r/2 rotations or all diagonal reflections are given 
nonzero probability ||^. Note that in the bulk, in step 1, we can restrict k to belong to 
{1,...,A^ — 1}, since pivot moves with pivot ljq are symmetry transformations. However, 
this is not the case in the presence of a surface, and it is thus important to include k = 0. 

A further requirement has to be satisfied in order to take into account the excluded 
region: 

4. If uj' ^ nff^uo), i.e. if S (uj'j) < for some i = 1, . . . , N, the move is rejected and the 
old walk is counted again. 

Because of the presence of the excluded region, ergodicity is not always satisfied. If S is 
a {d— l)-dimensional hyperplane, one can prove that axes reflections and 7r/2 rotations are 
enough to ensure ergodicity. However, even if ergodicity is satisfied, one may be worried by 
the fact that in order to explore the relevant phase space, the algorithm has to go through 
highly improbable configurations. For instance, if we denote by z the direction perpendic- 
ular to the hyperplane S, we expect the algorithm to be quite slow in updating the links 
at the beginning of the walk that are directed in the z direction. Indeed, transformations 
that have the pivot near the surface and that modify the z direction — moves involving the 
inversion of the 2;- axis, reflections with respect to the diagonals, 7i/2 rotations in the x = 
or y = planes — are unlikely to be accepted because of the excluded region. 

Therefore, it is important to study the dynamic behaviour of the algorithm as a function 
of the pivot point in order to determine whether the algorithm is still efficient. For a generic 
move based on a lattice symmetry transformation g applied in the pivot point we can 
define a local acceptance fraction fg^piv{N,k). In the scaling limit — > 00, /c — > 00 with 
k/N finite, the local acceptance fraction is expected to satisfy the scaling Ansatz 

fgMN, k) ^ N-P^--hg,p,, (k/N) , (3) 

where Pg,piv is the exponent which asymptotically governs the decay of the global acceptance 
fraction for N ^ 00 and hg^piy {k/N) is a scaling function. For small values of a = k/N 
the behaviour of the scaling function hg^pi^{a.) depends drastically on whether or not the 
transformation g preserves the z coordinate of the walk. 

In order of analyse quantitatively the dynamic behaviour of the algorithm, we group the 
pivot moves into equivalence classes. Two pivot moves which are based on the symmetries g 
and g' are said to be equivalent if there exists a global symmetry a — i.e. a transformation of 
the whole walk — which preserves the geometry of the system (i.e. a lattice symmetry which 
does not modify the z coordinates) such that, for every walk uj G ^'^(co'o) and for every pivot 
point ojk-i the walk which is obtained by applying a pivot move based on g' to the globally 
cr-transformed walk (y{(jj) in the pivot point a;^ = (7(0;^) and then transforming back via 
coincides with the walk obtained from uo by means of a pivot move based on the symmetry 
g with pivot point Uk- In formulae, if we denote by Gu the operator associated with the 
pivot move based on g applied on monomer fc, then g and g' are equivalent, if, for every fc. 
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Gk = cr^^ o G'l. o cr. For instance, consider the transformations g : (x, y, z) — >■ (x, z, y) and 
^f' : {x,y,z) — >■ (x, — — 2/). It is a simple matter to show that they are equivalent: It is 
enough to consider the 7i/2 rotation a : (x, y, z) — > (x, —z, y). 

In 3 dimensions the original 47 different moves are classified in 15 different classes. 
Moves belonging to a given class have exactly the same critical behaviour. Therefore, in 
the following we will study Eq. for different classes of moves. 

In Fig. |I| we report the acceptance fraction for the 15 classes as a function of a, for 

= 100. For small a, that is for pivot points that are near the surface, there are essentially 
three types of behaviour. Moves that involve the inversion of the z-axis have a very low 
acceptance rate: for a = 0.1 the acceptance of, say, reflections with respect to the 2;-axis 
is 0.04 and it drops further to 0.0045 for a = 0.05. Transformations that involve 7r/2 
rotations and diagonal reflections have a better behaviour although the acceptance still 
drops as q; — s> 0. Transformations that do not modify the z-direction are instead unaffected 
by the presence of the surface. 

Although the behaviour of the acceptance fraction as a function of a is radically different 
between SAW in the bulk and SAW in a half-space, the critical behaviour of the global 
acceptance fraction averaged over k is very similar. 

Denoting with fi^piv{N, k) the acceptance fraction of the class of moves i applied in the 
pivot point Uk, and with /^^™ its average over k, we expect /^^™ to vanish as the length 
of the walk increases as 

where the index i indicates a symmetry class, Aj is the corresponding leading correction- 
to-scaling exponent, and Ai, Bi are non-universal amplitudes. 

Before giving a quantitative estimate of the exponents Pi^piv, we try to give heuristically 
a rough estimate of the relative order of magnitude of the acceptance exponents in the bulk 
and in the presence of the surface. The argument which we use was already introduced in 
Ref. p] for SAW's in the bulk and we simply extend it to the case of SAW's in a half-space. 

If the two parts of the walk u [0, k] = {uq, . . . , Uk} and u [k, N] = {ujk, ■ ■ ■ , ^^n} are 
considered as independent, the local acceptance fraction for any applied lattice symmetry 
in the bulk would be given by 

AMN,k)^^^, (5) 

CkCN-k 

where is the number of walks of length k. For k — ^ oo, 

Ck ~ /fc^-' , (6) 
where fi is the critical fugacity and 7 ~ 1.1575(5) |T^, so that we have 

~ N'-^ , (7) 

and pi^pi^ ^ 7 - 1 ^ 0.16. 
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In the presence of the surface the same argument gives 



kMN,k)^—^^^^i-- (8) 

where Cfc(r) is the number of walks of length starting from r. For k oo, Ck{0) ~ fi'^k^'^^, 
where is a fugacity which is identical to that appearing in Eq. (^) and 7 a new critical 
exponent. In three dimensions, 7 ~ 0.68 |[T6|-[T8|. Assuming uJk to be at a (macroscopically) 



finite distance from the surface we can take CN-ki^k) oc CAr_fc, where c^^k scales according 
to Eq. §). Then, Eq. (D gives again f'P'" ~ A^^-^, i.e. pi^pi^ ^ 7 - 1 ^ 0.16. Thus, 
heuristically, we expect the acceptance exponents Pi^piv to be of the same order in the 
presence of the surface and in the bulk. This is confirmed by the numerical estimates. 

In order to compute the acceptance exponents Pi we considered walks of length = 
100, 200, 400, 800, 4000, 8000, 16000, 32000 and performed fits of the global acceptance 
fraction for each class of moves to a simple power law. In order to study the effect of 
the subleading terms appearing in (^) we performed different fits using only data for N > 
Nmin for increasing values of Nj^in = 100, . . . , 16000. We obtained as a result the effective 
exponents which are reported in table ^ and which approach pi as the lower cutoff A^min 
increases. 

Except for a few cases, which in the following we indicate with a star, the results are 
stable for walks of length A^ > 4000. Our best estimates for the different classes of moves 
are: 

1. (a) z-axis inversion: ppi^j = 0.0919(47)*, 
b) X or y-axis inversion: ppiy = 0.0925(11); 

a) ±tt/2 rotation in yz [or zx] planes: Ppiy = 0.10202(81), 

b) ib7r/2 rotation in the xy plane: Ppiv = 0.10002(60); 

a) vr rotation in yz [or zx] planes: Ppi^ = 0.1338(12), 

b) vr rotation in the xy plane: Ppi^ = 0.1329(17); 

a) diagonal reflection in the yz [or zx] planes: ppi^ = 0.08875(94)*, 

b) diagonal reflection in the xy plane: Ppiti = 0.0887(14)*; 

a) diagonal reflection in the yz [or zx] plane and x [resp. y] axis rcfl.: ppiv — 

0.13024(83), 

b) diagonal reflection in the xy plane and z-axis inversion : ppi^ = 0.1327(23); 

a) ±tt/2 rotation in the yz [or zx] planes and x [resp. y] axis refl.: ppi^ = 0.13556(76), 

b) ib7r/2 rotation in the xy plane and z-axis inversion: ppiy = 0.13371(69); 

7. 3-axes reflection: ppiy = 0.1637(17); 

8. diagonal reflection in the yz or zx planes and ib7r/2 rotation in the xy plane: ppi^j = 
0.11811(77); 

9. two diagonal reflections, one in the yz or in the zx plane and the other in the xy plane: 
Ppiy = 0.11813(94). 
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Here we have grouped the equivalence classes in the presence of the surface in 9 groups. 
These groups correspond to the 9 equivalence classes for the algorithm in the bulk, where 
one can consider global transformations a that do not preserve the z coordinate. From 
the results we observe that the exponent Ppj^, depends only on the equivalence class in 
the hulk. In agreement with the above-reported heuristic argument, the presence of the 
surface has no influence on the acceptance exponents of the different classes of moves, but 
only on the shape of the scaling functions hipiy{a) introduced in Eq. On the other 
hand, one should observe that the acceptance exponents which correspond to different bulk 
equivalence classes are different, with the only exception of classes 8. and 9., which can be 
hardly distinguished even in the presence of the surface. One should also notice that the 
exponent is larger for 3-axes reflections than for diagonal reflections. As it has been already 
noticed in Ref. this can be understood intuitively. Indeed, two subwalks on the opposite 
sides of a pivot point tend to be directed and occupy on average opposite regions of space. 
We can imagine, for instance, that the two subwalks occupy two opposite octants which 
touch at the pivot point. The subwalk uj [k, N] is not moved in the opposite octant by pivot 
moves based, for instance, on a one-axis reflection, a diagonal reflection or a 7r/2 rotation, 
while reflection will move the whole subwalk in the opposite octant. Therefore, it 

should have a higher probability of rejection. 

The stability of most of the fits with A^min seems to rule out the possibility that all 
equivalence classes have the same exponent and that the discrepancies are due to residual 
corrections to scaling. Different symmetries have apparently different acceptance expo- 
nents. If this is correct, the global acceptance exponent ged over all transforma- 
tions would coincide with the exponent of the pivot move which, in the limit N ^ (yo, has 
the highest probability of being accepted, that is Ppi^ ~ 0.089. This exponent is lower than 
that reported in Ref. which was obtained by averaging over all transformations. We 
mention that if we also perform the group average we obtain ppiy ^ 0.116(1). The group 
average exponent is only slightly higher than what can be found using the data reported 
in Ref. ||^ for the pivot algorithm in the bulk, for which p ^ 0.1 13Q. 

Now we can use our best estimates of the acceptance exponents for testing our scaling 
Ansatz (|^) and determining the behaviour of the scaling functions hi^piy{a) for a = k/N —>■ 

0. We will show that for a ^ 

hi,pUa) ~ a^' , (9) 

with Qi < for the "good" moves which do not change the z coordinate of the walk, and 
qi > otherwise. 

In Fig. we report fi^piv{N, k) ■ {N/lOOY'-p'" = hi^piy{a) for 4 different classes of moves: 

1. (b), 3.(b), 9. and 3. (a). They are representative of the different types of behaviour for 
a — ^ 0. One observes very good scaling, the data for different A^ falling one on top of 
each other, except at small values of a for those moves for which hi^piy is nonzero at small 

^ We have been informed that the acceptance fractions / reported in Ref. ||] for d = 3 (but not for 
d = 2) are in error because the program performed with probabihty 1 /48 an identity move (which is always 
accepted). Therefore, data there should be corrected by the map / (48/47) (/— 1/48). This mistake was 
guessed by Tom Kennedy on the basis of his own pivot simulations (see Ref. jl^] ) and was confirmed by 
Madras and Sokal's examination of their program. Aware of this problem, we recomputed the acceptance 
fraction exponent and found p ~ 0.113. 
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a in the range of values we have considered. Indeed, for moves l.(b) and 3.(b) the 
estimates increase as increases: for = 100 we would estimate /ii.(b),pjt)(0) ~ 0.92 and 
^3.(b),pw(0) ~ 0.84, while for = 32000 we have /ii.(b),pii,(0) ~ 1.12 and h3,(b),piv{0) ^ 1.10. 
For the moves of class 9., instead, the curve for A^ = 100 is slightly higher than that for 
A^ = 4000, 32000: for A^ = 100 we would estimate /i9.,pi^(0) ^ 0.15, while for A^ = 4000, 
32000 we would obtain /?.9.,pi^(0) ~ 0.07, although the data for A^ = 32000 seem to be even 
lower. Clearly, for small a, there are significant corrections to scaling and indeed we will 
now show that /ii.(b),piD(a) and /i3.(b),pi?;(a) diverge for a —>■ 0, while /i9.,pj^(0) = 0. 

Consider first moves in classes l.(b) and 3.(b) and in general all the moves preserving the 
z-coordinate of the walks. They are only marginally affected by the presence of the surface. 
The only effect is that, since the surface induces a monomer depletion near the surface, 
a move applied in a pivot point ujk with small k has higher probability of being accepted 
than the same move applied in UN-k, explaining the slight asymmetry of the scaling curves 
hi^piv{ot). Moreover, we expect the local acceptance fraction fi,piv{N, k) to remain finite as 
A^ increases with k fixed and small. By comparing with the scaling Ansatz (j^), it follows 
that hi^piy{a) ~ a'^^-p^", i.e. qi = —pi^piv and the scaling function diverges for a ^ 0. 
This is confirmed by the numerical results: indeed, /ii.(b),pTO(«) and /i3.(b),piD(«) increase as 
N ^oo. 

Let us now consider the classes of moves which change, but do not invert, the z- 
coordinate, for instance class 9. In order to understand the behaviour of hi^piv{a) for a — > 
we must compute the local acceptance probability for k fixed and small. For this purpose, 
we have considered the acceptance fraction averaged over the pivot points k < fcmax, i-e. 

kMN, A;,nax) = T -T /9-.P-(^' ^) " (^0) 

fcmax + i 

For fixed fcmax, it decays faster than the global acceptance fraction. For instance, for 
^max = 20, for the moves that involve n/2 rotations or diagonal reflections in the {y,z), 
{z,x) planes (classes 8. and 9.), we have numerically found /8.,9.„piD(A^, 20) ~ ]\f-0A9^ 
will be shown in the following that this exponent is close to the one which characterizes 
the probability that a generic bulk A^-step SAW grafted at the surface belongs to Qf^^uo). 
Since we expect that /8.,9.,p™(A^, ^max) ~ /8.,9.,pto(^) ^max) for small k, the scaling function 
h8.,9.,piv{c() is expected to vanish with gg-.g. ~ 0.49 — P8.,9.,piv ~ 0.37. Similar behaviour is 
expected for all classes of moves ( 2. (a), 4. (a), 5. (a), 6. (a), 8. and 9.) which do not invert 
the z axis. 

The behavior of moves like l.(a) that invert the z-axis is much worse. For instance, we 
performed a long run for A^ = 800, in which such moves were attempted on the first 20 
monomers approximately 21000 times, and none of the attempts was successful. Fitting 
directly /ii.(a)(a) for a = k/N < 0.1 for A^ = 100,800,4000, one finds gi.(a) = 3.8(3). 

It can be useful to get heuristically a rough idea of the order of magnitude of the 
exponent qi,[a). Let us denote with mN{uj,k) the maximum elongation in the z direction 
of the monomers following the pivot point 

mjq{(jj,k)= max u^, (11) 

i=k+l,...,N 
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where a;f is the 2;-coordinate of the monomer cUj. If we denote with uj^ the ^-coordinate of 
the pivot point, the proposed move does not give rise to colhsions with the surface if the 
inequahty 

> (12) 

is satisfied. Thus, if the pivot point is one of the first vertices following the grafted end, 
Eq. (|T2|) states that the move may be successful only for those walks that have nearly all 
their monomers within few lattice planes from the surface and therefore we expect that 
h{a) vanishes for a ^ 0. It is known (see Ref. [^]) that for oo, z oo and ( = z/N" 
fixed and small, the probability distribution Pn{z) of the fraction of monomers which lay 
at a distance z from the surface scales as 

Pn{z)^^X''"- (13) 

Therefore, the fraction W{z) of monomers that lay within a distance z from the surface 
scales as 

W{z)= [ Pn{z') dz' ^ (14) 







If k is small but already in the scaling regime, we expect that ~ k'^ and therefore 
W{2uj^) ~ k/N = a. Therefore, since the average fraction of monomers laying in the strip 
z < 2ujI is a, the average fraction of walks such that all their monomers lay in the strip is 
less than a. It follows that q\,(a) > 1. 

Clearly, the pivot move is inefficient when the pivot is near the surface. Therefore, 
we expect the dynamics of observables that strongly depend on the behavior of the walk 
near the surface — for instance, the number of monomers on the surface — to be much slower 
than that of global observables. Also, note that Tg^p becomes rapidly large, making it very 
difficult to thermalize the system using only the pivot algorithm. In the following Section 
we discuss a second nonlocal move that solves the problems we have discussed. 



3 The Cut- and- permute move 

In this Section we want to consider a different nonlocal move that is able to modify the 
walk near the surface. By adding it to the pivot algorithm we will obtain an algorithm 
without the quasi-ergodicity problems we have discussed above. Of course, we do not want 
to destroy the optimal dynamic behaviour of the pivot algorithm, and thus we want to 
introduce a move for which the mean CPU time per succesful move scales simply with A^. 
The move we introduce here consists in cutting the walk in two parts, and in rebuilding 
the walk in reverse order. 

The elementary cut-and-permute move works as follows: 

1. Given an A^-step SAW uo G VLf^{ujQ) , choose with probability Pc a cut-point ujc with 
c G {1, 2, ... , A^— 1}. The point ujc divides the walk in two subwalks uj^ = {ui, . . . , Uc} 
and u"^ = {uc, ■ ■ ■ , ujn}- 
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2. Cut the walk in the cut-point and rebuild it arranging the two subwalks and uj"^ in 
reverse order. The resulting walk u' has vertices at positions 




"^0 + ^c+i ~ 



for < 2 < - c, 
ioiN-c+l<i<N. 



(15) 



3. The proposed move is accepted if the resulting walk uj' is self- avoiding and does not 
intersect the excluded region. 

It is easy to see that the move satisfies detailed balance as long as Pc = Vn-c- 

It is also possible to incorporate in the cut-and-permute move the "good" pivot trans- 
formations, i.e. those that do not change the 2;-coordinate. We thus define an "improved" 
cut-and-permute move as follows (the move is illustrated in Fig. |]): 

1. Given an A^-step SAW uo G Vt^ , choose with probability pc a cut-point ujc with 
c G {1,2,..., A^— 1}. The point divides the walk in two subwalks = {ui, . . . , Uc} 
and u"^ = {uJc+i, . . .,ujn}- 

2. Choose with probability Pg an element g belonging to the symmetry group of the 
plane which is parallel to the surface. In three dimensions the symmetry group is 
the dihedral group D4 in the z = plane, whose 8 elements are ±7r/2 rotations, vr 
rotation, axis inversions, diagonal reflections, and the identity. 

3. Apply the chosen symmetry to the subwalk u'^ to obtain the subwalk u'^' = {uc, uJc + 

gi.^c+1 - Wc), • • • , + gi^^N - ^^c)}- 

4. Cut the walk in the cut-point and rebuild it arranging the two subwalks u'^ and cj^' 
in reverse order. The resulting walk uj' has vertices at positions 



5. The proposed move is accepted if the resulting walk uj' is self-avoiding and does not 
intersect the excluded region. 

The improved cut-and-permute move does not satisfy detailed balance. Indeed, the walk uj' 
obtained from uj with a cut-and-permute move in the cut-point uJc cannot be transformed 
back to the original walk uj by any move applied in cutv-c, since the subwalk of uj following 
the cut-point has changed its orientation in space and its orientation cannot be modified in 
the second cut-and-permute move. However, the move still leaves the probability measure 
invariant as long as pc = Pn-c- Indeed, if we perform a move based on g in the point 
Uc and then a move based again on g in the point ciJtv-c; we obtain the walk uj" = g [uj) 
corresponding to the application of symmetry g to the original walk. If we denote by 
p{uj uj') the probability that the move applied to uj gives as result uj', from what we have 
described it follows that p{uj — > uj') = p{uj' giy))- Thus, 




for < 2 < A^ - c, 

for A^ - c + 1 < i < A^. 



(16) 



1 



(17) 



UJ 



UJ 
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as required. In the following we will consider the improved cut-and-permute move, but, as 
we shall show explicitly, nothing would change by using the simpler version. 

As it has been done for the pivot algorithm, we can define equivalence classes of cut-and- 
permute moves and study the acceptance fraction fi^cp{N, c) as a function of the position 
c of the cut-point and its average /^^^^ over all points c. For large N, we expect a scaling 
behavior of the form 

/^^f ~ N-P^". (18) 

We have not added an index i to Pcp because, as we shall see, this exponent does not depend 
on the equivalence class. 

To give an estimate of the exponent pcp, we can use again a heuristic argument in which 
we consider the two subwalks in which the original walk uj is divided at the cut-point ujc 
as independent. We also assume that the probabihty that the attempted move is accepted 
is the product of the probability -Pi(c) that the concatenation of the two subwalks gives a 
walk which does not have self-collisions times the probability P2{c) that the subwalk lj2, 
transformed under a D4 symmetry and translated with its first vertex on the surface does 
not intersect the forbidden region. Repeating the argument presented for the pivot case, 
we have 

Pi{c) = J^, (19) 

CcCN-c 

which, by averaging over c, gives Pi ~ ]\[-i'y-^) jy-o. 1575(5) _ rjj^g second probability is 
given by 

P2ic) = ^ , (20) 

CN-c 

which, after averaging, gives P2 ~ N~^'^~'^^ ~ j\f-0A8_ r^^ie global acceptance fraction of 
the cut-and-permute move is expected to scale as fcp = Pi ■ P2 ~ N'^^p ~ jY-27+1+7 ^ 
jY-0.63 'jj-^jg argument is likely to give an overestimate of the exponent pcp for the following 
reasons. The original subwalk ui which is grafted at the surface is directed, so that Pi(c) 
is underestimated. Moreover, the walk uj2 is expected to have a residual directionality in 
the +z direction so that the probability that 002 does not intersect the forbidden region is 
expected to be larger than P2{c). Table contains numerical estimates of the acceptances 
/^'* for different classes of moves in D4 at different values of A^, averaged over the cutting 
point c. We have fitted the data for > A'min to a power law considering increasing 
values of A^min- Systematic errors due to corrections to scaling should become negligible for 
Amin 00. The values of fits with A^min > 4000 indicate that the fits are stable, with 
the exception of two cases. These two cases are denoted with a star in the following list 
containing our best estimates for the exponent pcp'- 

1. identity Pcp = 0.4969(35); 

2. diagonal reflections Pcp = 0.4903(29); 

3. ±tt/2 rotations Pcp = 0.4911(21); 

4. vr rotation pcp = 0.4873(40) '^; 

5. 1-axis reflection pcp = 0.4929(13)*. 
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At variance with the pivot case the results for the different equivalence classes are compat- 
ible within error bars: the transformation g seems to play little role. If we average over all 
symmetries we obtain pcp = 0.4922(20). Note that, as expected, this result is somewhat 
lower than the heuristic estimate reported above. 

The exponent pcp is significantly higher than ppi^ so that cut-and-permute moves will 
be more rarely accepted. Numerically, we find that a cut-and-permute move is accepted 
every ~ 2 (resp. ~ 20) succesful pivot moves for N = 100 (resp. N = 32000). How- 
ever, they represent the moves that most effectively change the conformation of the walk 
near the grafted end and thus they will play an important role in decorrelating the walk 
configurations. 

As for the pivot case, it is interesting to study the local acceptance fraction as a function 
of c and A^. For N ^ oo, ( = c/N fixed we expect a scaling form 

/,p(Ar,c)^Ar-P-/i,,(C), (21) 

where hcp{C) is a scaling function which encodes the dependence on the different cut-points. 
In Figure § we report the function hcp{C) = fcp{N, c) ■ (A^/100)^''*' for 5 values of A^ using the 
exponent pcp given by our fits. There is no difference in the local behaviour of the 5 classes 
of equivalent moves, so the function we report refers to the total acceptance obtained by 
averaging over the equivalence classes. 

Note the very good scaling: the points fall on top of each other for all values of A^. 

Let us analyze now the scaling behaviour of the acceptance at fixed c. We define the 
observable 

^ J2fcp{N,c)+ Yl ' (22) 

c=l c=N- 



fcp{N, 



2c„ 



and study its scaling with A^ as A^ — > oo. For Cmax = 20 a power-law fit gives fcp{N, 20) ~ 
jy-o.oi 'jj^jg jg jyg^ g^j^ effective exponent and indeed it decreases in modulus as the mini- 
mum value of A^ considered in the fit is increased. The results are therefore compatible with 
fcp{N, 20) ~ const for A'min > 4000. Such a result is easy to understand. Indeed, suppose 
first that A^ 

— Cmax ^ c < A'^. Then the move is accepted if co''[0, c'], c' — A'^ — c, does not 
intersect the surface — which happens with a probability independent of A^ — and does not 
intersect uj'[c', N], which, for large A^, should be roughly independent of A^. Thus, fcp{N, c) 
should be constant for A^ large. If 1 < c < Cmax, note that fcp{N,c) = fcp{N,N — c), 
to conclude the argument. In conclusion, the cut-and-permute move has a good scaling 
behavior exactly in the region in which the pivot moves behave badly. The addition of 
the cut-and-permute move in the hybrid dynamics is therefore crucial in speeding up the 
dynamics of the first steps of the walk. 

As in the pivot case, since the local acceptance fraction converges to a constant for 
small or large fixed values of (, the scaling function hcp{C) increases as ('P^p for — > 
(resp. (1 — Q'P^p for ( —>■ 1.) This suggests that in step 1. one can take pc different from 
only for ( near and 1, i.e. set pc = if aA^ < c < A^(l — a). For instance, by taking 
a = 1/10, we would approximately increase the acceptance by a factor 1.7. 
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4 Computational complexity 



In order to study the effective dynamic behaviour of the algorithm in CPU time, it is 
important to determine the scahng behaviour of the CPU time needed to generate pivot 
and cut-and-permute successful moves. 

The CPU time depends on the data structure that is used. For the walk we consider 
a sequentially allocated linear list and in order to check for self-collisions and collisions 
with the surface we use a hash table as described in Ref. [Q. This allows to insert a single 
monomer in an average CPU time of order one. In the following we will define a CPU time 
unit as the CPU time which is needed to upgrade a single monomer position. 

In order to optimize the algorithm, it is important to mimimize the CPU time spent in 
a failed attempt. This requires a careful choice of the order in which the new walk uj' is 
rebuilt. 

For the pivot algorithm, Ref. [§] suggested to build the new walk and perform the self- 
avoidance check in the following order: cu^, t^^+i, • • •jCo'fc+i, cUfc-i, • • • till all monomers 
have been checked. The motivation was that the probability of self-intersections is higher 
for those monomers which are close to the pivot point. However, in the presence of a 
surface a move that changes the z-coordinate can also fail because the new walk intersects 
the surface. 

Intersections with the surface occur with higher probability for those monomers of the 
subwalk uj[k,N] which are more likely to undergo a large displacement, while monomers 
belonging to the grafted subwalk uj [0, k] are not modified and satisfy the geometrical con- 
straint automatically. In order to keep into account this effect, for the moves that change 
the ^-coordinates, we can use a different strategy: we insert the monomers of the proposed 
walk ui' in the hash table in the following order: tu^, c<j^, uj'^^j^^, c<j^_i, ^'n-ii ■ ■ •> ^'k+v ^k-v 
ujN-i, • • • , till every monomer has been checked once. As a measure of the CPU time spent 
by the algorithm we can use the number of walk monomers I{u}') that are checked in the 
move. If the move is successful, clearly I{uj') = N, while for a failed attempt 

I{uj') = mm{i : uj' [k - i,k - 1] n u' [k, k + i] $ 

OTuj'[N-i,N] n S ^(lSoTuj'[k-i,k + i] n S^0}, (23) 

where, following Ref. [§, we denote the subwalk {<ax(i,o)' <ax(i,o)+i' " " " <in(i,7V)}' ^i^h 
the symbol u' [i,j] and the half-space z < with S. We want now to evaluate the average 
time spent in checking a failed move that is expected to scale as 

Tr'^""'^'^ (/(^'))lfailed~Ar^-^ (24) 

where the average over all failed moves is taken. An estimate of ypiy can be obtained by a 
heuristic argument. Inserting monomers in the hash table in the way we described before, 
the probability that I{uj') exceeds the value i can be estimated from the scaling of the 
acceptance fraction. Indeed, we expect 

Prob{/(co'') > i} = Prob{co''[A; — i, k + i]is SAW and does not intersect the surface, 

and uj'[N — i,N] does not intersect the surface } ~ i~Ppi^ . (25) 
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Therefore, the average time spent in checking a failed move is expected to scale like 

N 

j.piv Jailed ^ Prob{/ > z} ~ N^-P^^- , (26) 

i=0 

so that Upi^ = 1 — ppi^. Estimates of j'P^'^^f"-^'-^'^ are reported in Table ^ for the different 
classes of pivot moves and different values of A^. We performed fits of the form (^4]) in 
order to determine ypi^, including each time only data with > A^min (in Table ^, A^min is 
the value A^ reported in the first column). We observed that in all cases i/pi^ increases with 
A^min and thus our results are probably a lower estimate of Upiy. Here, we report the results 
with A^min = 4000 and indicate with a star those cases in which the is still too large: 

1. z-axis inversion: Upi^ = 0.8750(16), 

X or y-axis inversion: ypj^, = 0.8234(17)*; 

2. ib7r/2 rotation in yz [or zx\ planes: Hpiv — 

0.8554(13)*, 

ib7r/2 rotation in the xy plane: ypiv = 0.8205(23); 

3. vr rotation in yz [or zx] planes: ypiv = 0.82756(97)*, 
7T rotation in the xy plane: ypiv — 0.7989(46); 

4. diagonal reflection in the yz [or zx] planes: ypiv = 0.8635(22), 
diagonal reflection in the xy plane: ypiv = 0.8470(25)*; 

5. diagonal reflection in the yz [or zx] plane and x [resp. y] axis refl.: ypiv = 0.8381(13)*, 
diagonal reflection in the xy plane and z-axis inversion : ypiv = 0.8342(16); 

6. ib7r/2 rotation in the yz [or zx] planes and x [resp. y] axis refl.: ypiv = 0.8250(18), 
ib7r/2 rotation in the xy plane and z-axis inversion: ypiv = 0.8339(23); 

7. 3-axis reflection: ypiv = 0.8138(35); 

8. diagonal reflection in the yz or zx planes and ib7r/2 rotation in the xy plane: ypiv = 
0.8476(20)*; 

9. two diagonal reflections, one in the yz or in the zx plane and the other in the xy plane: 
ypi„ = 0.85121(67)*. 

If we average over aU equivalence classes we obtain i/piv = 0.8428(13). These results should 
not be trusted too much, and in many cases the correct estimate is probably higher. How- 
ever, the important thing that emerges is that i/piv ^ 1 — Ppiv, a result that will be used 
below. 

We also computed 7"^™'/"*'^'^ for the moves that change the z-coordinate when we insert 
the points as proposed in Ref. i.e. without modifying the order in which we insert 
monomers into the hash table because of the presence of the surface. We found that, by 
using this method, the average number of points inserted in the hash table before detecting 
a failure is approximately 1.4 times larger. The difference is more significant for pivot 
points Uk with k < N/2 — for k/N the difference is a factor of two — while for k > 0.6 
the two methods are equivalent. However, this has no influence on the exponent ypiv, which 
is approximately the same in both cases. 
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At this point we can estimate the average CPU time spent in a pivot move 

rp /-| rpiv\rppiv, failed . rpivrppiv,succ 

where we have taken into account that ypi^ < 1 —ppiv Notice also that more time is spent 
in a successful move than in an unsuccesful one. 

Let us now consider the cut-and-permute move. As in the previous case we tried to 
build the new walk in the most efficient way in order to reduce the CPU time necessary to 
detect a failure. Unlike the pivot case, for the cut-and-permute move the heuristic argument 
indicates that the dominant cause of failure is given by intersections with the surface. As 
it is illustrated in Fig. ^, the cut-point Uc divides the walk in two subwalks, ui [0, c], which 
is grafted at the surface, and u [c, N] . The move transforms u [c, N] via a lattice symmetry 
g which preserves z and translates the transformed subwalk grafting its first vertex at the 
surface. The subwalk u[0, c] is instead translated in order to join the free end of the grafted 
subwalk. The only possibility of hitting the surface comes from the new grafted subwalk 
u' [0, N — c] and the monomers which are more likely to intersect the forbidden region are 
the ones following the new grafted vertex. Collision between monomers are instead more 
likely to occur at the joining of the two subwalks. 

As for the pivot move, in order to check at the same time the two regions where failures 
are more likely to occur, we inserted the monomers of the proposed new walk uj' in the 
hash table in the following order: u;^.^, uJq, cj^v-c-d ^n-c+v " " ^i-i' ^N-c-i^ ^N-c+i^ 

In this way the CPU time needed to check a proposed walk is 

I{uj') = mm{i ■.uj'[N-c-i,N-c-l]nLu'[N-c+l,N-c + t]^$ 

OTUj'[0,i] n S ^dloTuj'lN -c~i,N -c + i] n S ^ 0}, (28) 

if u' is rejected and I{uj') = if it is accepted. The average time spent in detecting a 
failed move is expected to scale as 

^e^Jarled^ (/(o.') ) Ifailed ~ ^T^^ • (29) 

We can estimate i/cp heuristically as before. With the above procedure the probability of 
using CPU time which is longer than i goes as 

Prob{/(c<;') > i} = 

Prob{ct;'[A^ — k — i, N — k + i]is SAW and u!'[0, i] does not intersect the surface } 
~ i-P^^ , (30) 

and the average total time spent in a failure scales as 

N 

r^epjailed ^ ^ p^o6{/ > i} ~ N^-^^^ , (31) 
i=0 

SO that Ucp = 1 —Pcp- The data for j'^'f"-^^^'^ in the cut-and-permute move and the estimates 
of the i/cp exponent obtained from data with < Nmin, with Nmin appearing in the ffist 
column, are reported in Table From those data we have computed our final estimates as 
we did in the pivot case obtaining: 



Clearly 
(27) 
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1. diagonal reflections yep = 0.4870(19); 

2. identity yep = 0.4788(70); 

3. ±7r/2 rotations yep = 0.4921(13); 

4. TT rotation yep = 0.4933(43)*; 

5. 1 axis reflection y^p = 0.4869(13)*. 

The group average yep — 0.4890(15) is slightly lower than 1 — ~ 0.5078(20) as it is also 
happened for the pivot move. 

In order to confirm the role played by intersections with the surface, we also measured 
/(a;') when monomers are inserted in the hash table in the order uj'j^_e, i-^N-e+i^ '^n-c-i^ 
■ ■ without inserting at the same time the monomers near the grafted end. In this case 
we find yep ~ 0.9 which is very similar to the exponent for the pivot move. The proposed 
procedure is thus very inefficient. 

Finally, we can estimate the average CPU time spent in a cut-and-permute move. 
Clearly 

rp /-i fCp\ rpCp, failed , rCp rpCp,SUCC 

-i-N,cp — JnI'-'-n '^Jn'-'-n 

~ 1 • N^-^ + N-P-^ ■ N ~ ATi-Pcp , (32) 

where we have taken into account |/cp ^ 1 — Pep- Note also that TN,ep/TN,piv ~ N^^p'^p^" ~ 
N"'^-^, so that the average CPU time for successful move scales identically. 

5 The hybrid cut-and-permute algorithm: autocorre- 
lation times 

In this Section we want finally to define our hybrid algorithm that consists in performing 
both pivot and cut-and-permute moves. More precisely, the algorithm is specified by a 
number < g < 1. The hybrid algorithm works by performing a pivot move with prob- 
ability q and a cut-and-permute move with probability (1 — q). Since the pivot alone is 
ergodic, the full algorithm is ergodic. 

In order to understand the dynamic behavior we will distinguish three different classes 
of observables: (a) global observables that depend on the global shape of the SAW, for 
instance, the end-to-end distance or the radius of gyration; (b) global surface observables 
that depend on the interaction of the SAW with the surface, for example, the number of 
monomers lying on S; (c) local observables that depend on local properties of the SAW. 

Let us analyse at first the behaviour of global observables. Both pivot and cut-and- 
permute are nonlocal moves and can decorrelate global observables in few accepted moves, 
but the cut-and-permute move cannot modify a class of observables, for instance the z- 
coordinate of the free endpoint. We denote with Oz the class of observables which are 
updated only in the pivot dynamics and with Ox,y all the others. For a polymer of length 
N the autocorrelation time Tint,Oz is therefore of order l/(g • fN,piv), while Tint,o^^y is of 
order l/fN,max, where fN,max = max{5 • fN,piv, (1 - • fN,cp}- Since in the limit N ^ oo 
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the acceptance fraction of the cut-and-permute move decreases faster than the acceptance 
fraction of the pivot move, if q is constant, the scahng behaviour of the autocorrelation 
time Tint,giob of any global observable is determined by the pivot dynamics only, i.e. 

TinLglob ~ Tintfi. ~ T-i„i,o.,, ~ " A^^*"'" , (33) 

SO that Zgiob ~ Ppiv ~ 0.10. Table ^ reports the autocorrelation times for different values 
of for the z component of the free endpoint obtained setting q = 1/2. The dynamic 
exponent z, which is obtained by fitting the data to a power law, is z ~ 0.15, in reasonable 
agreement with the prediction. 

Global surface observables will be changed in a few successful cut-and-permute moves 
and thus we expect 

r.„M-/ ~ (1 - ?)^''^^ (34) 

so that Zsurf ~ Pep ~ 0.49. Unfortunately, we have not measured surface quantities in our 
simulation and therefore we are not able to test this prediction. Nonetheless, we believe 
that the correct value of z is not too far from pcp. 

Finally, let us consider a local observable, for instance the average number of monomers 
at which the SAW makes a 90° turn. Such a quantity is of order A^ with a variance of order 
\/N. By using a standard random-walk argument we obtain 



loc ~ mm{q-^N^+P^^% (1 - qy^N^+P^^} ~ N^+Pv^^ ~ N^-^\ (35) 



since Pcp < Ppiv Since local moves should be the slowest ones we expect Texp ^ A^^'^. 

Finally, we estimate the behaviour of the autocorrelation time for global observables 
in CPU units. Since the average CPU time spent in performing one move of the cut-and- 
permute algorithm is proportional to 

T = q-Tpi^ + {l-q)-T^p^q- N^-Pp^^ + (l _ g) . N^-p^p ^ g . n^-Pv.. ^ (35) 
the autocorrelation time of global observables in CPU units scales as 



r. 



where, in the absence of a numerical estimate, we have used our heuristic estimate, Zsurf ~ 
Pcp- The behavior of global observables is nearly optimal, while for surface observables 
the behavior is not as good. If one is interested in computing surface observables, it 
is possible to improve the algorithm by increasing the frequency of the cut-and-permute 
moves. Indeed, if we scale q~^ ~ NP'^^'P''''' — essentially we keep the relative frequency of 
successful pivot and cut-and-permute moves constant — then our previous arguments give 

rint,glob ~ A^P^f , Tint,surf ~ A^^^^ (39) 

The time per move scales now as N^'P'^p^ so that 

^CPU ^CPU Ar /An\ 



int,glob 
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In this case, the dynamics is optimal for both types of observables. 

Note that one cannot obtain the same result using the pure pivot algorithm. Indeed, the 
autocorrelation time for global surface observables would be proportional to the number 
of attempts occurring between two accepted pivot moves changing the ^-coordinate of the 
walk and applied in a pivot point uok with k small. Because of the monomer depiction near 
the surface, indeed, one does not expect that a move at a pivot point far from the grafted 
end would affect macroscopically the number of monomers near the surface. In Sec. 2 we 
have shown that among the z-changing moves only the ones which do not invert z have 
a small probability of being accepted for small k. This probabihty decays as with 
X 0.49. Therefore, if we choose the pivot point with uniform probability, we get 

and the autocorrelation time in CPU units scales as 

^SSr - N^:^^'-^^^^ ^ N'-' . (42) 

This indicates clearely that the pivot algorithm alone is inefficient for simulating polymers 
grafted at an inpenetrable interface. 
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Figure 1: Acceptance fraction for the different classes of pivot moves as a function of the 
rescaled variable a = k/N, for = 100. Three kinds of behaviour are visible at small 
values of a. At a = the acceptance fraction vanishes for 5 classes of moves transforming 
z — > —z, it is small for 6 classes corresponding to n/2 rotations and diagonal reflections in 
the {y, z) and {z, x) planes, while it is not affected by the presence of the surface for the 
remaining 4 classes of moves which do not modify the z coordinate. The correspondence 
between symbols and classes of lattice symmetries as defined in Sec. 2 is the following: 
l.(a) 2;-axis inversion (□), l.(b) x or y-axis inversion (*), 2. (a) ±7r/2 rotation in yz [or zx\ 
planes (■), 2.(b) ±7r/2 rotation in the xy plane (V), 3. (a) tt rotation in yz [or zx] planes 
(+), 3.(b) TT rotation in the xy plane (x), 4. (a) diagonal reflection in the yz [or zx] planes 
(o), 4.(b) diagonal reflection in the xy plane (a), 5. (a) diagonal reflection in the yz [or zx] 
plane and x [resp. y] axis refl. (•), 5.(b) diagonal reflection in the xy plane and 2;-axis 
inversion (▼), 6. (a) ±7r/2 rotation in the yz [or zx] planes and x [resp. y] axis rcfl.(A), 
6.(b) ±7r/2 rotation in the xy plane and 2;-axis inversion (0), 7. 3-axes reflection (♦), 8. 
diagonal reflection in the yz or zx planes and ±7r/2 rotation in the xy plane (pentagons), 
9. two diagonal reflections, one in the yz or in the zx plane and the other in the xy plane 
(full pentagons). 
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Figure 2: Scaling functions in the variable k/N for the pivot acceptance fraction. The 4 
classes of moves considered are, in order of increasing global acceptance: 3. (a), 9., 3.(b) 
and l.(b). 
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Figure 3: A cut-and-permute move: here c— 17 and the symmetry g is an inversion of the 
X-axis. 
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u 


32000. 


26075 ± 00128 








l.(b) 


{x,y,z) 


{-x,y,z)\/ {x,- 


y,z) 




100. 


0.73933 ±0.00010 


0.08937 ±0.00006 


240.326 


6 


200. 


0.69716 ±0.00018 


0.09035 ± 0.00009 


38.684 


5 


400. 


0.65530 ±0.00039 


0.09102 ± 0.00017 


17.865 


4 


800. 


0.61704 ±0.00031 


0.09167 ±0.00024 


2.573 


3 


4000. 


0.53241 ±0.00048 


0.09251 ±0.00113 


1.988 


2 


8000. 


0.49970 ±0.00012 


0.09511 ±0.00216 


0.002 


1 


16000. 


0.46785 ± 0.00092 


0.09542 ± 0.00712 


0.000 





32000. 


0.43791 ±0.00198 








2. (a) 


{x, y, z) (-z, y, x) V (z, y, -x) V (a;, - 


z,y) V {x,z,-y) 


100. 


0.49426 ±0.00013 


0.10258 ±0.00008 


92.048 


6 


200. 


0.46014 ±0.00107 


0.10203 ± 0.00010 


2.764 


5 


400. 


0.42768 ± 0.00005 


0.10202 ±0.00010 


1.642 


4 


800. 


0.39837 ±0.00014 


0.10191 ±0.00019 


1.142 


3 


4000. 


0.33805 ±0.00021 


0.10202 ± 0.00081 


1.121 


2 


8000. 


0.31508 ±0.00010 


0.10337 ±0.00161 


0.188 


1 


16000. 


0.29391 ± 0.00146 


0.10670 ±0.00786 


0.000 





32000. 


0.27296 ± 0.00062 
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2.(b) 


{x,y,z) 


{y,-x,z) V {-y,x 






100. 


0.67272 ± 0.00006 


0.09676 ±0.00005 


1269.317 


6 


200. 


0.62958 ± 0.00065 


0.10012 ±0.00011 


28.757 


5 


400. 


0.59014 ±0.00005 


0.10018 ±0.00011 


6.794 


4 


800. 


0.55070 ±0.00006 


0.10043 ±0.00015 


1.166 


3 


4000. 


0.46879 ± 0.00077 


0.10002 ±0.00060 


0.663 


2 


8000. 


0.43679 ± 0.00025 


0.09978 ± 0.00067 


0.010 


1 


16000. 


0.40768 ± 0.00070 


0.10003 ± 0.00270 


0.000 





32000. 


0.38037 ± 0.00028 








3. (a.) 


[.r.y.z) ^ 


(•'■■-.'/■--) V (-.i:./y 


■---) 




100. 


0.34586 ±0.00017 


0.13240 ±0.00015 


138.561 


6 


200. 


0.31323 ±0.00014 


0.13106 ±0.00019 


7.738 


5 


400. 


0.28620 ± 0.00078 


0.13156 ±0.00034 


4.618 


4 


800. 


0.26141 ±0.00013 


0.13158 ±0.00035 


4.535 


3 


4000. 


0.21167 ±0.00010 


0.13383 ±0.00117 


0.457 


2 


8000. 


0.19300 ±0.00021 


0.13462 ±0.00269 


0.350 


1 


16000. 


0.17556 ±0.00049 


0.13054 ±0.00740 


0.000 





32000. 


0.16037 ± 0.00069 








3.(1)) 


(:r. 


,j.z)^i-.v.-y.z} 






100. 


0.58549 ± 0.00055 


0.13128 ±0.00016 


37.268 


6 


200. 


0.53739 ± 0.00008 


0.13145 ±0.00016 


6.018 


5 


400. 


0.49089 ± 0.00020 


0.13147 ±0.00036 


6.014 


4 


800. 


0.44774 ±0.00010 


0.13086 ±0.00046 


1.565 


3 


4000. 


0.36295 ± 0.00032 


0.13291 ±0.00172 


0.044 


2 


8000. 


0.33084 ± 0.00094 


0.13242 ±0.00358 


0.020 


1 


16000. 


0.30205 ±0.00158 


0.13367 ±0.00964 


0.000 





32000. 


0.27532 ±0.00114 
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4. (a) 


{x,y,z) {x,z 


, y) V {x, -z, -y) V {z, y, x) V (- 


-z,y,-x) 


100. 


0.50437 ±0.00017 


0.09358 ±0.00010 


165.230 


6 


200. 


0.47158 ±0.00072 


0.09278 ±0.00012 


40.909 


5 


400. 


0.44132 ± 0.00008 


0.09277 ±0.00012 


39.705 


4 


800. 


0.41398 ±0.00001 


0.09291 ±0.00014 


34.737 


3 


4000. 


0.35583 ± 0.00027 


0.09094 ±0.00049 


16.959 


2 


8000. 


0.33377 ±0.00039 


0.08875 ± 0.00094 


9.538 


1 


16000. 


0.31270 ± 0.00022 


0.08538 ±0.00144 


0.000 





32000. 


0.29473 ± 0.00021 








4.(1)) 


(,r,/y.c)- (/y..r,-jV(-.(y,-,r,,:) 


100. 


0.68479 ± 0.00026 


0.09004 ±0.00012 


154.206 


6 


200. 


0.64473 ± 0.00098 


0.09132 ±0.00016 


13.171 


5 


400. 


0.60772 ± 0.00023 


0.09140 ±0.00016 


6.351 


4 


800. 


0.57025 ±0.00015 


0.09133 ±0.00019 


5.922 


3 


4000. 


0.49255 ± 0.00035 


0.09107 ±0.00081 


5.808 


2 


8000. 


0.46179 ±0.00022 


0.08866 ±0.00141 


1.507 


1 


16000. 


0.43448 ± 0.00041 


0.09672 ± 0.00671 


0.000 





32000. 


0.40631 ± 0.00185 








5, (a.) 


(:r. z) ^ (-.v. z.y) V (-.,;. - . -ij) V (, 


.. -IJ..V) V (--, -.(/, -,r) 


100. 


0.40080 ± 0.00015 


0.13372 ±0.00013 


709.276 


6 


200. 


0.36448 ± 0.00073 


0.13174 ±0.00015 


21.659 


5 


400. 


0.32989 ±0.00004 


0.13170 ±0.00015 


4.219 


4 


800. 


0.30091 ± 0.00033 


0.13124 ±0.00039 


2.576 


3 


4000. 


0.24312 ±0.00029 


0.13024 ±0.00083 


0.710 


2 


8000. 


0.22245 ± 0.00042 


0.13062 ±0.00179 


0.654 


1 


16000. 


0.20294 ±0.00013 


0.12917 ±0.00253 


0.000 





32000. 


0.18556 ±0.00030 
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5.(b) 


{x, y, z) ^ (y, X, -z) V {-y, -x, -z) 


100. 


0.32210 ± 0.00024 


0.13234 ±0.00030 


46.991 


6 


200. 


0.29528 ± 0.00052 


0.13115 ±0.00045 


33.813 


5 


400. 


0.26675 ± 0.00019 


0.13021 ±0.00048 


4.807 


4 


800. 


0.24413 ±0.00015 


0.13098 ±0.00063 


1.138 


3 


4000. 


0.19781 ±0.00022 


0.13268 ±0.00228 


0.532 


2 


8000. 


0.18073 ±0.00048 


0.13566 ±0.00469 


0.006 


1 


16000. 


0.16458 ±0.00098 


0.13657 ±0.01264 


0.000 





32000. 


0.14971 ± 0.00096 








6. (a) 


(x.fj.z) ^ {-.V.Z.- 


-/yjV(-,r,-,-.,/y)V(- 


---,-,';,.':) 


V (.--, -.y, -:r) 


100. 


0.39576 ± 0.00008 


0.13914 ±0.00009 


149.867 


6 


200. 


0.35778 ± 0.00084 


0.13683 ±0.00021 


5.301 


5 


400. 


0.32505 ± 0.00030 


0.13683 ±0.00021 


5.299 


4 


800. 


0.29599 ±0.00002 


0.13687 ±0.00022 


3.607 


3 


4000. 


0.23732 ±0.00012 


0.13556 ±0.00076 


0.284 


2 


8000. 


0.21601 ± 0.00022 


0.13525 ±0.00149 


0.226 


1 


16000. 


0.19652 ± 0.00035 


0.13363 ±0.00372 


0.000 





32000. 


0.17913 ± 0.00034 








6.(1)) 




z)^{-tj.A:.-i)\J{ji.-.r.-z) 




100. 


0.31919 ± 0.00020 


0.13629 ±0.00011 


108.973 


6 


200. 


0.28836 ±0.00025 


0.13629 ±0.00012 


108.963 


5 


400. 


0.26430 ± 0.00003 


0.13642 ±0.00012 


35.570 


4 


800. 


0.23944 ± 0.00023 


0.13422 ±0.00040 


2.428 


3 


4000. 


0.19272 ±0.00010 


0.13371 ±0.00069 


1.603 


2 


8000. 


0.17576 ±0.00007 


0.13506 ±0.00127 


0.002 


1 


16000. 


0.16007 ± 0.00034 


0.13522 ±0.00406 


0.000 





32000. 


0.14575 ±0.00027 
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7. 



{x,y,z) {-x,-y,-z) 



100. 

200. 

400. 

800. 
4000. 
8000. 
16000. 
32000. 



0.27382 ± 

0.24357 ± 
0.21512 ± 
0.19260 ± 
0.14818 ± 
0.13222 ± 
0.11810 ± 
0.10580 ± 



0.00011 

0.00012 
0.00000 
0.00041 
0.00015 
0.00014 
0.00063 
0.00067 



0.16839 ±0, 

0.16438 ±0, 
0.16223 ±0, 
0.16332 ±0, 
0.16370 ±0, 
0.16129 ± 0, 
0.15868 ±0, 



00019 

00025 
00027 
00095 
00174 
00404 
01196 



1141.917 


6 


516.155 


5 


1.987 


4 


0.554 


3 


0.488 


2 


0.054 


1 


0.000 






(x, y, z) {y, z, -x) V {y, -z, x) V {-y, z, x) V {-y, -z, -x) 
{x, y, z) {z, X, -y) V {z, -x, y) V {-z, x, y) V {-z, -x, -y) 



100. 

200. 

400. 

800. 
4000. 
8000. 
16000. 
32000. 



0.42766 ± 
0.39272 ± 
0.36094 ± 
0.33180 ± 
0.27392 ± 
0.25247 ± 
0.23265 ± 
0.21389 ± 



0.00006 
0.00011 
0.00019 
0.00008 
0.00016 
0.00009 
0.00043 
0.00052 



0.12071 ±0.00007 


198.212 


6 


0.11972 ±0.00011 


66.816 


5 


0.11894 ±0.00015 


11.374 


4 


0.11872 ±0.00017 


1.545 


3 


0.11811 ±0.00077 


0.871 


2 


0.11911 ±0.00151 


0.282 


1 


0.12131 ±0.00442 


0.000 






9. (.T, ?y, z) (y, z, x) V {y, -z, ~x) V (-y, -z, x) V {-y, z, -x) 

(x, y, z) (z, x, y) V (z, —x, —y) V (— z, —x, y) V (— z, x, —y) 



100. 

200. 

400. 

800. 
4000. 
8000. 
16000. 
32000. 



0.42794 ± 
0.39344 ± 
0.36165 ± 
0.33281 ± 
0.27512 ± 
0.25351 ± 
0.23351 ± 
0.21520 ± 



0.00003 
0.00081 
0.00006 
0.00010 
0.00026 
0.00017 
0.00065 
0.00040 



0.12050 ±0, 
0.11872 ±0, 
0.11871 ±0, 
0.11821 ±0, 
0.11813±0, 
0.11821 ±0, 
0.11781 ±0, 



00008 
00017 
00017 
00026 
00094 
00138 
00485 



138.509 
8.153 
7.110 
0.022 
0.013 
0.007 
0.000 



Table 1: Acceptance fraction and acceptance exponent for different equivalence classes of 
pivot moves. 
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N 


In "JN 


Pep ± Spcp 




DF 


1. 




100. 


0.24511 ±0.00008 


0.44881 ± 0.00021 


3544.599 


6 


200. 


0.18075 ±0.00071 


0.47947 ± 0.00056 


81.376 


5 


400. 


0.13252 ±0.00001 


0.47979 ± 0.00057 


49.763 


4 


800. 


0.09563 ±0.00012 


0.48524 ±0.00105 


12.206 


3 


4000. 


0.04397 ±0.00007 


0.49686 ± 0.00352 


0.218 


2 


8000. 


0.03111 ±0.00015 


0.49507 ± 0.00751 


0.146 


1 


16000. 


0.02214 ±0.00019 


0.50161 ±0.01870 


0.000 





32000. 


0.01564 ±0.00015 








2. 


{x,y) {y,x) V {~y,-x) 


100. 


0.24517 ±0.00001 


0.46012 ± 0.00023 


2179.892 


6 


200. 


0.18130 ±0.00027 


0.47995 ±0.00051 


321.318 


5 


400. 


0.13259 ±0.00015 


0.48458 ± 0.00059 


65.902 


4 


800. 


0.09570 ±0.00007 


0.48753 ± 0.00069 


1.363 


3 


4000. 


0.04373 ± 0.00012 


0.49032 ±0.00295 


0.415 


2 


8000. 


0.03116 ±0.00006 


0.49210 ±0.00467 


0.176 


1 


16000. 


0.02219 ±0.00011 


0.49719 ±0.01301 


0.000 





32000. 


0.01572 ±0.00012 








3. 


ix,y) V iy,-x) 


100. 


0.24518 ±0.00006 


0.46855 ±0.00011 


16481.760 


6 


200. 


0.18142 ±0.00039 


0.48350 ± 0.00016 


753.217 


5 


400. 


0.13251 ±0.00003 


0.48369 ± 0.00016 


633.232 


4 


800. 


0.09576 ±0.00004 


0.48822 ±0.00024 


2.753 


3 


4000. 


0.04386 ±0.00015 


0.49106 ±0.00211 


0.915 


2 


8000. 


0.03112 ±0.00001 


0.49005 ± 0.00236 


0.003 


1 


16000. 


0.02217 ±0.00021 


0.49086 ±0.01458 


0.000 





32000. 


0.01577 ±0.00005 









4. 


{x,y) {-x,-y) 


100. 


0.24507 ±0.00016 


0.46150 ±0.00035 


2088.088 


6 


200. 


0.18193 ±0.00047 


0.48081 ± 0.00056 


153.010 


5 


400. 


0.13278 ±0.00007 


0.48182 ±0.00058 


93.048 


4 


800. 


0.09579 ± 0.00007 


0.48739 ± 0.00083 


3.736 


3 


4000. 


0.04365 ±0.00015 


0.48733 ±0.00402 


3.736 


2 


8000. 


0.03121 ±0.00007 


0.49257 ±0.00698 


2.894 


1 


16000. 


0.02234 ±0.00014 


0.52741 ± 0.02164 


0.000 





32000. 


0.01550 ±0.00021 








5. 


{x,y) {-x,y) V {x,-y) 


100. 


0.24513 ±0.00012 


0.46950 ±0.00015 


10482.119 


6 


200. 


0.18120 ±0.00055 


0.48043 ± 0.00019 


2075.983 


5 


400. 


0.13249 ±0.00003 


0.48056 ±0.00019 


2003.417 


4 


800. 


0.09554 ± 0.00000 


0.48857 ±0.00026 


21.472 


3 


4000. 


0.04384 ±0.00009 


0.49292 ±0.00136 


10.893 


2 


8000. 


0.03117 ±0.00007 


0.49062 ±0.00231 


9.382 


1 


16000. 


0.02207 ±0.00002 


0.48175 ± 0.00371 


0.000 





32000. 


0.01580 ± 0.00004 









Table 2: Acceptance fraction and acceptance exponents for equivalent classes of cut-and- 
permute moves. 
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N 


rppiv^ failed _|_ ^rppiv , j ailed 


Upiv ^ ^JJpiv 




DF 


l.(a) 


{x,y 


,z) {x,y,-z) 






100. 


12.48970 ±0.01355 


0.85395 ± 0.00038 


1572.322 


6 


200. 


21.88020 ±0.06886 


0.86694 ± 0.00051 


116.838 


5 


400. 


38.71450 ±0.04355 


0.86837 ±0.00053 


22.278 


4 


800. 


70.38660 ±0.13222 


0.87022 ± 0.00079 


12.277 


3 


4000. 


283.17500 ±0.72023 


0.87498 ± 0.00163 


1.072 


2 


8000. 


516.57200 ±2.50767 


0.87807 ±0.00381 


0.264 


1 


16000. 


954.17300 ±8.43733 


0.87159 ±0.01318 


0.000 





32000. 


1745.83000 ± 3.99927 








l.(b) 


{x,y,z) 


{-x,y,z)\/ {x,-y,z 


) 




100. 


16.83500 ±0.00583 


0.78863 ± 0.00030 


818.622 


6 


200. 


28.39370 ±0.14783 


0.80608 ±0.00077 


210.977 


5 


400. 


47.96640 ±0.17382 


0.80994 ±0.00087 


122.163 


4 


800. 


82.41200 ±0.31495 


0.81591 ±0.00116 


61.053 


3 


4000. 


299.93200 ±0.67181 


0.82335 ±0.00166 


21.853 


2 


8000. 


526.85700 ± 1.76380 


0.83244 ±0.00315 


10.320 


1 


16000. 


926.25200 ±2.84146 


0.84895 ± 0.00603 


0.000 





32000. 


1668.36000 ±4.73391 








2. (a) 


(x, y, z) {-z, y, x) V {z, y, -x) V {x, -z, y) V (x, z, -y) 


100. 


13.90580 ±0.00651 


0.82973 ± 0.00021 


2599.558 


6 


200. 


24.36420 ±0.07489 


0.84812 ±0.00042 


112.260 


5 


400. 


43.11250 ±0.11340 


0.84883 ± 0.00043 


45.278 


4 


800. 


77.02470 ± 0.01481 


0.84900 ± 0.00043 


35.698 


3 


4000. 


300.49700 ±0.49023 


0.85542 ±0.00127 


6.812 


2 


8000. 


540.37900 ± 1.11283 


0.85976 ±0.00212 


0.329 


1 


16000. 


976.35600 ± 7.32634 


0.86609 ±0.01124 


0.000 





32000. 


1779.62000 ±3.74866 
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2.(b) 


{x,y,z) 


{y,-x,z) V {-y,x,z) 




100. 


14.49390 ±0.01393 


0.78652 ± 0.00038 


637.741 


6 


200. 


24.40480 ± 0.06407 


0.79844 ±0.00068 


193.153 


5 


400. 


41.18000 ±0.12866 


0.80573 ± 0.00092 


55.520 


4 


800. 


71.24760 ±0.20412 


0.80964 ±0.00125 


33.978 


3 


4000. 


258.58000 ±0.49761 


0.82051 ± 0.00229 


2.083 


2 


8000. 


455.17500 ±0.91544 


0.82495 ± 0.00396 


0.198 


1 


16000. 


807.80900 ± 3.65702 


0.82059 ±0.01058 


0.000 





32000. 


1426.69000 ± 8.22899 








3. (a.) 




i;.-(y,--) V {-x.tj.- 






100. 


11.52150 ± 0.00774 


0.81006 ±0.00012 


18940.759 


6 


200. 


19.36200 ±0.05615 


0.81986 ±0.00014 


2917.879 


5 


400. 


33.04310 ±0.03296 


0.82019 ±0.00014 


2558.777 


4 


800. 


57.17680 ±0.01268 


0.82104 ±0.00015 


2056.303 


3 


4000. 


212.82400 ±0.03815 


0.83872 ± 0.00044 


207.907 


2 


8000. 


382.59400 ±0.16969 


0.82756 ±0.00097 


42.738 


1 


16000. 


674.83200 ±0.71975 


0.84356 ± 0.00263 


0.000 





32000. 


1210.97000 ± 1.79389 








3.(1,) 


(:r,/y. 


z) ^ (-x.-!j.z) 






100. 


13.49420 ± 0.01353 


0.74416 ± 0.00047 


4523.560 


6 


200. 


21.84060 ±0.15401 


0.77124 ±0.00062 


111.280 


5 


400. 


35.90370 ±0.00864 


0.77146 ±0.00062 


82.659 


4 


800. 


60.44850 ± 0.12562 


0.77967 ±0.00123 


22.250 


3 


4000. 


209.65400 ±0.73649 


0.79887 ±0.00461 


3.622 


2 


8000. 


363.04100 ±0.81002 


0.81292 ±0.00905 


0.368 


1 


16000. 


636.07500 ±4.66641 


0.83143 ±0.03183 


0.000 





32000. 


1131.86000 ±23.55230 
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4. (a) 


{x, y, z) {x, z, y) V {x, -z, -y) V [z, y, x) V {-z, y, - 


x) 


100. 


13.55210 ± 0.00593 


0.83932 ±0.00032 


1054.187 


6 


200. 


23.77720 ±0.02946 


0.85113 ±0.00051 


163.158 


5 


400. 


42.22660 ±0.05254 


0.85557 ±0.00063 


26.317 


4 


800. 


75.99580 ±0.10340 


0.85770 ±0.00086 


12.732 


3 


4000. 


301.14100 ±0.61872 


0.86352 ±0.00219 


4.372 


2 


8000. 


544.65900 ± 1.36295 


0.86954 ± 0.00363 


0.029 


1 


16000. 


996.21800 ±6.31179 


0.86773 ±0.01119 


0.000 





32000. 


1817.89000 ±8.13467 








4.(1)) 


[.r.y.z) 


[y- ■>■.:-) V (-(/.-.I 


,-) 




100. 


13.92810 ±0.01710 


0.79980 ±0.00034 


2770.326 


6 


200. 


23.31840 ±0.07178 


0.81499 ± 0.00046 


370.328 


5 


400. 


40.02870 ±0.10436 


0.81686 ± 0.00048 


194.954 


4 


800. 


69.26540 ± 0.04385 


0.81769 ±0.00050 


145.693 


3 


4000. 


256.45000 ±0.83364 


0.83938 ± 0.00209 


30.983 


2 


8000. 


451.01700 ±0.58539 


0.84708 ± 0.00251 


0.871 


1 


16000. 


809.00600 ±2.64035 


0.85312 ± 0.00694 


0.000 





32000. 


1461.39000 ± 5.17075 








5, (a.) 


[x.y.z) ^ (-,i:,-,/y) V ( 


-/y)V(-.-/y 


.,r)V (-,-., - 


-.r) 


100. 


12.19490 ± 0.00138 


0.79610 ± 0.00009 


13776.870 


6 


200. 


20.54390 ±0.04020 


0.81876 ±0.00022 


977.370 


5 


400. 


35.10850 ±0.04365 


0.82039 ± 0.00023 


407.213 


4 


800. 


60.83820 ±0.02967 


0.82179 ±0.00025 


180.301 


3 


4000. 


225.88200 ±0.43779 


0.83565 ±0.00109 


10.462 


2 


8000. 


400.64600 ± 0.27076 


0.83814 ±0.00134 


0.074 


1 


16000. 


716.28200 ±0.46864 


0.83497 ±0.01175 


0.000 





32000. 


1277.72000 ± 10.36998 
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5.(b) 


ix,y,z) - 


^ (2/,a;,-2;)V(-?/, - 


x,-z) 




100. 


10.63500 ±0.00391 


0.79204 ±0.00013 


6828.502 


6 


200. 


17.72970 ±0.02629 


0.81120 ±0.00027 


621.314 


5 


400. 


30.14010 ±0.02065 


0.81429 ±0.00031 


179.047 


4 


800. 


52.36940 ±0.10904 


0.82540 ±0.00101 


46.110 


3 


4000. 


195.91700 ±0.06318 


0.83425 ±0.00165 


0.147 


2 


8000. 


349.47400 ± 0.61497 


0.83322 ± 0.00358 


0.042 


1 


16000. 


621.76000 ±4.43418 


0.83566 ±0.01243 


0.000 





32000. 


1109.64000 ±5.36154 








6. (a) 


(:r,/y.-.j ^ {-.v.z.-ij) 


V(-,r,-c,/y)V(-c, 


-/y,.,;)V(,-.. 


-u- -■'■) 


100. 


12.31430 ±0.01259 


0.79843 ± 0.00032 


3167.891 


6 


200. 


20.63680 ±0.05582 


0.80658 ± 0.00035 


228.953 


5 


400. 


35.42490 ±0.00217 


0.80680 ± 0.00035 


180.458 


4 


800. 


61.23020 ±0.08125 


0.81473 ±0.00075 


38.007 


3 


4000. 


225.32100 ±0.29390 


0.82502 ±0.00184 


0.293 


2 


8000. 


399.46000 ± 0.66259 


0.82396 ±0.00337 


0.152 


1 


16000. 


705.99200 ±3.17945 


0.82742 ± 0.00948 


0.000 





32000. 


1252.79000 ± 5.99259 








6.(1)) 


z) -r (-.(y,.i;. - -) V (.(y, - 


i:,-,-.) 




100. 


10.72070 ±0.01194 


0.79288 ±0.00041 


3124.323 


6 


200. 


17.77960 ±0.03092 


0.81265 ±0.00057 


677.940 


5 


400. 


30.26860 ±0.12081 


0.82221 ± 0.00069 


50.629 


4 


800. 


52.47530 ±0.04430 


0.82291 ± 0.00070 


25.697 


3 


4000. 


195.23300 ±0.83783 


0.83388 ± 0.00233 


1.349 


2 


8000. 


346.10300 ±0.87541 


0.83597 ±0.00295 


0.001 


1 


16000. 


617.91800 ±4.27681 


0.83572 ±0.01101 


0.000 





32000. 


1102.83000 ±3.55154 
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7. 


{x,y,z) {-x,-y,-z) 


100. 


10.43374 ± 0.00562 


0.72192 ± 0.00026 


6189.621 


6 


200. 


16.87000 ± 0.15782 


0.79349 ± 0.00097 


262.240 


5 


400. 


27.59870 ± 0.08282 


0.79487 ± 0.00098 


164.803 


4 


800. 


46.24860 ± 0.00737 


0.79758 ±0.00100 


24.244 


3 


4000. 


163.67700 ± 1.03384 


0.81382 ± 0.00354 


1.338 


2 


8000. 


285.89400 ± 1.67814 


0.81769 ± 0.00538 


0.421 


1 


16000. 


501.00000 ± 3.85448 


0.82535 ±0.01297 


0.000 





32000. 


887.75500 ±4.12887 








8. 


{x, y, z) (y, z, -x) V (y, ~z, x) V {-y, z, x) V (-y, -z, 
{x, y, z) -> (z, x, ~y) V (z, ~x, y) V (-z, x, y) V (-z, -x, 


-x) 

-y) 


100. 


12.46370 ± 0.00513 


0.80822 ± 0.00021 


7715.820 


6 


200. 


21.30210 ±0.00692 


0.81994 ± 0.00026 


2238.419 


5 


400. 


37.03280 ±0.00907 


0.82865 ± 0.00033 


242.017 


4 


800. 


65. 16220 ± 0.05450 


0.83399 ± 0.00052 


66.482 


3 


4000. 


247.15100 ± 0.43455 


0.84409 ±0.00140 


5.681 


2 


8000. 


441. 53500 ±0.69148 


0.84757 ± 0.00205 


0.269 


1 


16000. 


792.69500 ± 3.48873 


0.85116 ± 0.00723 


0.000 





32000. 


1429.98000 ±3.41729 








9. 


[x, y, z) -> (y, z, x) V (y, -z, -x) V (-y, -z, x) V (-y, z, 
(x, y, z) "> (z, x, y) V (z, -x, -y) V (-z, -x, y) V (-z, 


-x) 

-y) 


100. 


12.44090 ± 0.00677 


0.82676 ± 0.00012 


18185.883 


6 


200. 


21.30550 ±0.01874 


0.83542 ± 0.00014 


2970.316 


5 


400. 


37.06200 ± 0.06948 


0.83729 ± 0.00014 


715.471 


4 


800. 


65.09910 ±0.01040 


0.83740 ± 0.00014 


629.169 


3 


4000. 


246.75200 ± 0.26749 


0.84924 ± 0.00050 


24.712 


2 


8000. 


442.15800 ± 0.32010 


0.85121 ± 0.00067 


4.541 


1 


16000. 


793.87800 ± 1.71463 


0.85794 ± 0.00323 


0.000 





32000. 


1438.86000 ±0.83247 









Table 3: Average time spent in a failed pivot move and scaling-behaviour exponent t/piy 
for the equivalence classes of pivot moves. The check of self-avoidance was performed as 
explained in section ^ 
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N 


rpcp^j atied I crpcp^j ailed 


yep ± Sycp 




DF 


1. 




100. 


iU. 14144 zt U.UUIOO 


U.4( z ( y ± U.UUUzi 


0(^Q 7Q/1 
Zoo. / o4 


D 


200. 


lo.youzu =t u.UioU / 


U.4( yo4 ± U.UUU4/ 


Of^ /I Q1 
Z0.4ol 


c 



400. 


1 o /I HQ /in -u n noHQi 

iy.4Uo4U zt U.UZUol 


U.4oU ( 4 ± U.UUUD4 


1 ni A 

lO.UlD 


A 

4 


800. 


Z ( .Uol4U =t U.UloZ / 


U.4oi4o ± U.UUU / D 


1 OQO 






4000. 


0O.0UU4U =t U.Uol i 1 


U.4yUzo ± U.UUozo 


4.D0D 




z 


8000. 


CO 1=^/1 1=^1 n -u n 07*^nfi 
oZ.040lU =t U.Z / Dyo 


U.4 ( ooU ± U.UU ( Uz 


1 007 

i.zy / 


1 
1 


16000. 


llO.lolUUzt l.UOU/o 


n /I ^700 4- n ni oc;7 
U.40 / yy ± u.uiyt) < 


n nnn 

U.UUU 


u 


32000. 


1 "iQ 00+1 fil 97fi 








2. 


y) 


i^y, X) V [ y, X) 






100. 


lU.141ol =t u.uuiy/ 


n /171 7Q 4- n nnni 
U.4( i ± U.UUU it) 


101/3 oon 
lyiD.zzu 


D 


200. 


1 Q nc:Qon -u n nni on 
lo.yooyU =t U.UUlzU 


u.4(oiy ± u.uuuio 


001 0^^7 

ZZl.ZD / 


cr 



400. 


in Q^iOQnj- n c\aka i 


U.4ooO( ± U.UUUOD 


7 QOQ 
i .6Z0 


A 

4 


800. 


Zo.yOooU =t U.UU / 04 


U.4o4Uo ± U.UUUO / 


0.000 


Q 



4000. 


^Q. ^^i^^in-un 1 00 Qf; 
OO.DIDIU zt U.lZZoO 


n /I CAOQ 4- n nni oo 
u.4ooyo ± u.uuiyz 


n fi!^7 




z 


8000. 


Qo c\^A Qn -u n 0/1 Q1 fi 

oZ. U04oU zt U.Z4olo 


n /I con/i 4- n nnQC/i 

U.4oyU4 zt U.UUoo4 


n /1 7n 
u.4 i u 


1 
1 


16000. 


11/1 0/1/1 nn -u n osc;/i 7 

114.y44UU zt U.Zo04/ 


n /10Q7Q 4- n nn7fiQ 
U.4yo lo± U.UU ( oo 


n nnn 

U.UUU 


u 


32000. 


161 85000 + 78166 








3. 


{x,y) 


(-y,a;) V (y, -x) 






100. 


10.13721 ±0.00407 


0.47885 ± 0.00011 


584.104 


6 


200. 


13.97810 ±0.00776 


0.48137 ± 0.00016 


149.121 


5 


400. 


19.40490 ±0.02226 


0.48388 ± 0.00031 


57.209 


4 


800. 


26.98040 ±0.02724 


0.48550 ±0.00044 


30.702 


3 


4000. 


58.54480 ±0.06015 


0.49209 ±0.00129 


1.102 


2 


8000. 


82.38670 ±0.02772 


0.49019 ± 0.00250 


0.322 


1 


16000. 


115.68200 ±0.20951 


0.49841 ±0.01471 


0.000 





32000. 


163.41900 ± 1.63933 









4. 


{x,y) {-x,-y) 


100. 


10.14555 ±0.00415 


0.48117 ± 0.00006 


2964.627 


6 


200. 


13.94730 ±0.02649 


0.48198 ± 0.00006 


36.688 


5 


400. 


19.34220 ±0.00095 


0.48199 ± 0.00006 


22.650 


4 


800. 


27.00100 ±0.00421 


0.48226 ± 0.00010 


12.053 


3 


4000. 


58.61070 ±0.12987 


0.48721 ± 0.00257 


8.353 


2 


8000. 


81.96580 ±0.01518 


0.49328 ±0.00432 


5.292 


1 


16000. 


113.95800 ±0.70179 


0.52218 ±0.01329 


0.000 





32000. 


163.65800 ± 1.12058 








5. 


{x,y) 


^ {~x,y) V {x,~y) 






100. 


10.14238 ±0.00359 


0.47655 ± 0.00021 


344.868 


6 


200. 


13.93060 ±0.02251 


0.48196 ± 0.00037 


24.212 


5 


400. 


19.39390 ±0.01713 


0.48242 ±0.00041 


17.624 


4 


800. 


26.97070 ±0.04447 


0.48468 ± 0.00073 


3.912 


3 


4000. 


58.65150 ±0.07087 


0.48689 ±0.00134 


0.034 


2 


8000. 


82.21440 ±0.12537 


0.48668 ± 0.00231 


0.022 


1 


16000. 


115.14900 ±0.34679 


0.48752 ± 0.00610 


0.000 





32000. 


161.44300 ±0.47969 









Table 4: Average time spent in checking a failed cut-and — permute move for different 
euivalence classes of moves. The check of the proposed configurations is performed as 
described in section | 
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N 




z ^5z 


x' 


DF 


4000 


10.347 ±0.040 


0.1506 ±0.0029 


1.564 


2 


8000 


11.420 ±0.037 


0.1543 ±0.0041 


0.007 


1 


16000 


12.703 ± 0.087 


0.155 ±0.012 


0.000 





32000 


14.146 ±0.068 









Table 5: Autocorrelation times and dynamic exponent Zgiob for the observable z' 



36 



